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. Abstract 

^ \ We investigate the coherence properties of an atomic beam evaporatively cooled in a magnetic guide, 

assuming thermal equilibrium in the quantum degenerate regime. The gas experiences two-dimensional, 
transverse Bose-Einstein condensation rather than a full three-dimensional condensation because of the 
very elongated geometry of the magnetic guide. First order and second order correlation functions of the 
atomic field are used to characterize the coherence properties of the gas along the axis of the guide. The 
coherence length of the gas is found to be much larger than the thermal de Broglie wavelength in the 
strongly quantum degenerate regime. Large intensity fluctuations present in the ideal Bose gas model 
are found to be strongly reduced by repulsive atomic interactions; this conclusion is obtained with a 
one-dimensional classical field approximation valid when the temperature of the gas is much higher than 
its chemical potential, fc,sT ^> 

O ; Pacs: 03.75.Fi, 42.50.-p 

The first experimental achievements of Bose-Einstein condensates in atomic vapors Q ^, [3|, |J] have 
opened promising perspectives for atom optics: condensates constitute indeed atomic waves sources of much 
better coherence properties than the usual 'thermal' sources like the standard magneto-optical trap. These 
CO ' coherence properties have already been demonstrated experimentally: interferences experiments between 
■ two condensates have been performed at MIT || and at JILA || , the first order correlation function of the 
\q atomic field has been measured in Munich @, and suppression of density fluctuations (that is fluctuations 
| in the intensity of the atomic field) has been revealed by a measurement of the mean-field energy j|] and of 
three-body losses ||. By inducing a coherent leak of atoms out of trapped condensates several groups have 
succeeded in creating pulsed or quasi-continuous 'atom-lasers' [jOJ. 

For future applications the already realized 'atom-lasers' may suffer from the handicap of a low mean 
£h ' nux °f atoms: the condensates were not experiencing any continuous loading of atoms, so that the coherent 
i-^j ■ output of atoms terminated once the ~ 10 6 atoms of the condensate were leaked out. As the repetition rate 
of the whole sequence is limited by the time required to form a condensate by evaporative cooling (on the 
^ ' order of seconds) the resulting mean flux of atoms is < 10 6 atoms/s. Several proposals have been made to 



■4— > 

CO 



a 



O 

o 



refill the condensates with atoms in a continuous way [11| but they have to our knowledge not been realized 
yet. Recently we proposed a different scheme, based on the evaporative cooling of an atomic beam [p^ |. 
The goal of the present article is to predict the essential features of our 'continuous atom-laser' proposal. 
■ The paper is organized as follows. In section |l| we summarize the calculations of |TI| performed to determine 
the required length of evaporative cooling of the beam to reach the quantum degenerate regime. In section [2] 
we discuss the coherence properties of the beam once quantum degeneracy has been obtained, assuming that 
the atomic interactions are negligible; we find that the beam has large intensity fluctuations incompatible 
with expected coherence properties of an atom-laser. In section || we propose a model of a one-dimensional 
interacting Bose gas and we construct a classical field approximation to this model in the high temperature 
limit. We solve the classical field approximation using the formal analogy between functional integrals and 
quantum propagators: we find that the interactions between particles can dramatically reduce the intensity 
fluctuations of the beam. We conclude in section f|. 
*E-mail: yvan.castin@lkb.ens.fr 
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1 Evaporative cooling of an atomic beam 



1.1 A continuous injection of atoms in a magnetic guide 

In the experimental scenario considered at the Ecole normale superieure the continuous source of atoms is 
provided by standard laser cooling and trapping techniques taking place inside a cell. Atoms in the cell are 
captured, cooled and trapped in the x — y plane using a two-dimensional magneto-optical trap, that is with 
laser beams in the x — y plane and a dipolar magnetic field B oc xe x — ye y , where e x , y are unit vectors along 
x, y axes. The atomic motion along z is controlled with a standard moving molasses technique: the two 
counter-propagating laser beams along z have different frequencies so that the atoms are cooled around a 
non-zero mean velocity vq. 

The already cold atomic beam emerging from this set-up is sent towards a magnetic guide of axis z. The 
magnetic guide is produced by the superposition of a uniform magnetic field along z and a dipolar magnetic 
field oc xe x — ye y . This provides a transverse confinement of the atoms being in the right Zeeman sublevels, 
the atoms in the wrong Zeeman sublevels being not trapped or even expelled. The magnetic guide provides 
a trapping potential in the x — y plane with a harmonic bottom that we write as 

U(x,y) = ^muj 2 ± (x 2 + y 2 ) (1) 

where m is the atomic mass and u>± the transverse oscillation frequency of the atoms. 

The following parameters of the injected atoms are expected to be realistic for Rb atoms [12]. The 



initial velocity dispersion of the atoms is Avq = 20 cm/s corresponding to an initial temperature of 400 fiK; 
such a high temperature is obtained after spatial compression of the cloud usually performed to increase 
the collision rate in preparation of evaporative cooling. We take an injection velocity vq = 3Avq larger than 
the velocity spread so that the incoming atoms form a beam. Assuming an injected flux of atoms of 3 x 10 9 
atoms/s and an oscillation frequency co± = 2ir x 1 kHz we find an initial on-axis thermal density of 8 x 10 11 
atoms/cm 3 . The initial phase space density is 7 x 10~ 7 <C 1. The initial collision rate of atoms is related 
to the on-axis density no and to the collisional cross section a by 7^ = (2/yfn) no cr A^o- For the s-wave 
collisional cross-section of rubidium (er = 7.6 x 10~ 16 m 2 ) we get 7^ ~ 100 s _1 which is much smaller than 

(o) 

= 0.02 < 1. (2) 

The atoms have therefore the time to perform a full transverse oscillation in the trapping potential before 
experiencing a collision. 

1.2 Modeling of evaporative cooling 

We assume that the atoms are subject to evaporative cooling in the magnetic guide, e.g. by application of 
a z-dependent radio-frequency flipping the atoms to untrapped or expelled Zeeman sublevels when they are 
too far from z axis. 

The dynamics of evaporative cooling can be described by the classical Boltzmann equation on the phase 
space atomic density f(f,p) as long as the phase space density remains small. The collision terms in the 
Boltzmann equation are simplified by the assumption (justified for rubidium) that atomic interactions take 
place in the s-wave only and have a constant (momentum independent) total cross-section a = 8ira^ d (a^d 
is the scattering length). The effect of evaporative cooling is modeled by setting f(r,p) to zero in the 
domain x 2 + y 2 > A 2 (z) for two-dimensional evaporation, or in the domain x 2 > A 2 (z) for one-dimensional 
evaporation, where A(z) is an adjustable cut in position space. 

We have first looked for an approximate analytical solution of the Boltzmann equation, adapting to our 



geometry the truncated Gaussian ansatz put forward in [13]. As argued in [12] it is then more rigorous to 



restrict to a one-dimensional evaporation scheme. The ansatz then takes the following form: 

f(r,p) = /o(z)e-(- + ^)/fe T W)e-^-^)) 2 /( 2 ^^W)y (eevap ( 2 ) - €x ) (3) 
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where e x , e y are the sum of kinetic energy and harmonic trapping potential energy along x and y respectively 
and Y is the Heaviside function. The ansatz assumes a local thermal equilibrium with temperature T{z). 
The temperature depends on z only not on x and y as, due to Eq.(||), the mean free path of the particles 
~ Au/7 co ii is much larger than the spatial transverse extension of the gas (fcsT/(mw^_)) 1//2 ~ Av/uj±: 
transversally the gas is in the so-called collisionless regime. For the same reason the truncation of / in 
position space is replaced by a truncation in energy space, with 



r cvap 



(z) = -muj 2 ± A 2 (z). 



(4) 



Knowing the energy of the particle along x allows to calculate the maximal excursion of the trajectory along 
x (as atoms have in general the time to perform a full harmonic oscillation before experiencing a collision); 
if this maximal excursion exceeds A(z) the particle is evaporated. 

There are a priori three unknown functions of z in the ansatz Eq.(|3|): (i) the normalization factor fo(z) 
or equivalently the linear density pn n (z), (h) the mean momentum p{z) of the gas along z, and (iii) the 
temperature of the gas T(z). 

By multiplying Boltzmann's equation by (i) unity, (ii) the momentum p z , and (iii) the kinetic energy 
along z, p 2 /(2m), and by integrating over x,y,p x ,p y ,p z one gets three hydrodynamic type equations for 
Pim(z), p(z) and ksT^z). These three equations contain the usual equations expressing conservation of 
probability of momentum and of energy, plus extra terms describing the loss of particles, the change of 
momentum and energy due to the evaporation. 

We have solved numerically the hydrodynamic type equations in steady state, assuming that the z- 



dependence of the parameter A(z) is adjusted to maintain a z-independent ratio r\ 
the specific set of parameters of 



c evap 



[z)/(k B T(z)). For 

j |l . 1| we have to gain seven orders of magnitude on the phase space density 
to reach quantum degeneracy. The smallest spatial length of evaporation required is obtained for r/ ~ 5 as 
shown in figure |l[ it is equal to 7600 do where do = y/x/ (2noa) is the mean free path at the entrance of the 
magnetic guide, that is ~ 11 meters for the considered parameters. After evaporative cooling along these 
11 meters the flux of particles has been reduced by a factor 90 and the temperature has been decreased by 
a factor 4000. 
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Figure 1: For one-dimensional evaporative cooling, gain in phase space density as a function of the position z in the 
magnetic guide. The position is expressed in units of the initial mean free path do — y/if / (2no<r) . These curves are 
the numerical solutions of the hydrodynamic type equations in the case where rj = e cvap (z) / (kBT(z)) is fixed. 
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We have also performed a numerical simulation of the full Boltzmann equation by a Monte Carlo algo- 
rithm using macro-atoms [14|. The resulting numerical calculations take several days on a workstation. We 
have first simulated the one-dimensional evaporation scheme. In this way we have confirmed the accuracy of 
the predictions based on the ansatz Eq.(||). We have also performed simulations for the full two-dimensional 
evaporation. We have found typically that the length required to reach quantum degeneracy is reduced by a 
factor three as compared to one-dimensional evaporation, with the same loss of two orders of magnitude on 
the flux. For the specific set of parameters of §[0] the required evaporation length reaches the experimentally 
reasonable value of 4 meters. 



2 Coherence properties in the ideal Bose gas model 

We now assume that evaporative cooling has allowed to reach the quantum degenerate regime at a certain 
distance from the entrance of the tube. We have not performed any kinetic study of the approach of quantum 
degeneracy, and we will here simply assume that the gas is at thermal equilibrium in the frame moving at 
the mean velocity of the gas. Such an assumption is reasonable if the temperature remains significantly 
larger than the quantum fiw±_ of transverse oscillation of the atoms of the guide. If the temperature was 
much smaller than Tiu±/kB the transverse degrees of freedom of the gas would be frozen in the ground state 
of the transverse harmonic oscillator; the gas would become a free one-dimensional Bose gas along z, that 
could not thermalize as colliding identical particles in one dimension simply exchange their momenta. 

In this section we consider the model of the ideal Bose gas. The effects of atomic interactions are 
discussed in the next section. 



2.1 Transverse Bose-Einstein condensation 



Let us enclose the ideal Bose gas in a fictitious box of size L along z, with periodic boundary conditions. 
Transversally the gas is confined by the harmonic potential Eq.(|]). The one-particle eigenstates of the 
system are then labeled by three integers, the non-negative integers l x ,ly labeling the eigenstates of the 
harmonic oscillator along x and y, and the integer l z labeling the momentum along z: 



hk z 



^1 



In the grand canonical ensemble the mean occupation number of the single particle state I 
given by 



n(l,n) 



exp 



P 



2m 



1 



(5) 



(6) 



For convenience we have included the transverse zero-point motion energy in the chemical potential, so that 
\i varies from — oo to 0. 

It turns out that in our trapping geometry there is no Bose-Einstein condensation in the thermodynamical 
limit defined as L, N — > +oo with a fixed linear density N/L, N being the mean number of particles. 

Let us consider indeed a fixed value of L and let us define the maximal mean number of particles 
that can be put in all states but the ground state of the trap at a fixed temperature. A Bose-Einstein 
condensate forms in the ground state of the trap when ./V exceeds iV max . As each n(l,n) (for I ^ 0) reaches 
its maximal accessible value for /i = 0, iV max is given by 



(7) 



A lower bound on A^ ax is obtained by restricting the sum to the states transversally in the ground state of 
the harmonic oscillator: 

-l 



N' > V 

- 1 * max — / . 



eXP (^)" 1 



(8) 
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When L is large enough so that fcgT 3> /i 2 / (2mL 2 ) the exponential in the denominator can be expanded to 
first order, leading to 



/ r \ 2 
7T / L 



K^>- 3 [- x ) (») 

where we have introduced the thermal de Broglie wavelength 

A= (27rmA; B r) 1 /2- (10) 

One then realizes that -/V max grows faster than iV in the thermodynamical limit L — > +00, with N/L 
fixed. This is connected to the known fact that there is no Bose-Einstein condensation in a one-dimensional 
homogeneous Bose gas in the thermodynamical limit. 



On the other hand, transverse Bose-Einstein condensation [15| is taking place in our system. Let us 
calculate indeed the maximal number of atoms that can be put for a fixed temperature in the transversally 
excited single particle states: 

iV max= E E n (^ = °)' ( n ) 

(^,^(0,0) h 

Replacing in the large L limit the sum over l z by an integral and expanding l/(exp(x) — 1) in powers of 
exp(— x) we obtain: 



J * max — \ 

A S >1 



E s ~ 1 ' 2 { [1 " ex P i-'P*>±T a - l} - J 2 C(5/2). (12) 



The last equality is correct in the limit ksT S> Huj±. The function ( is the Zeta function of Riemann, and 
£(5/2) ~ 1.341. If the linear density of the gas pn n = N/L exceeds the critical value 

p£ = <»/l4(M) 2 C (5/2, (13 ) 

the excess of density will accumulate in the transverse ground state of the trap. At linear densities much 

(c) 

higher than the gas becomes almost monomode transversally: this is an interesting feature for atom 
optics applications as we have achieved in this way a Heisenberg limited transverse focalization of the beam. 

We can also consider the maximal on-axis density of atoms in the excited transverse states. In this 
case it is more convenient to label the eigenstates of the two-dimensional harmonic oscillator by the angular 
momentum M along z and the radial quantum number l T . Then the wavefunctions on z axis, that is in 
x = y = 0, have a squared modulus equal to zero if M 7^ and equal to mu>±/(Trti) for M = 0. We recall 
that the states with M = have an energy above the zero-point energy given by 2l r hoj±. We then obtain 
the maximal on-axis density of atoms in excited transverse states: 



L>1 L v - • 



(14) 



With the same algebraic transformations as for Eq.(12) we obtain 



(c) rnu; ± 1/2 

/'axis 7r?lA 

s>l 



1 — exp(— 2sf3huj±) 



where £(3/2) ~ 2.612 and we have used kgT 3> huj±. Transverse condensation therefore takes place when 
the usual Einstein's condition pA 3 = £(3/2) is satisfied on the axis of the trap! 
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2.2 In the quantum degenerate regime 

We assume now that the gas is in the strongly degenerate regime, with a linear density p\\ n larger than the 
critical value Eq. (|i~3| ) . The linear density of atoms in the excited transverse states has reached its saturated 
value p{?\ This implies that is much smaller than hu>± (fi can be replaced by zero for the transversally 
excited states); as k B T 3> hu>± one has also 

\fj,\<k B T (16) 

so that the occupation number of the absolute trap ground state / = 0, though not of order N, is much 
larger than unity: 

<l = 0) = 1 » 1. (17) 

exp(— pfjL) — 1 



Let us calculate the linear density of atoms in the transverse ground state of the trap: 

l^n(0,0,Z 2 ;/x) (18) 

h 

j Y, {exp[/?(^ 2 2 /(2m) - M)] - l}" 1 . (19) 



p {0) 
Hin 



L 

lz 



As (3\p\ <C 1 and the sum is one-dimensional the main contribution to the sum comes from states with 
kinetic energies on the order of \fi\. Expanding the exponential in the Bose formula to first order we obtain 
a Lorentzian approximation for the occupation number as function of momentum: 



ti 2 k 2 z /(2m) + \n 

Replacing finally the sum by an integral we obtain 



n(0, 0, l z ; n) ~ t2; 9;/ f ^ ; — ■ (20) 



This allows to express the chemical potential as function of density, when combined with the relation 

Pto = &+&. (22) 

Note that such a calculation would fail for a two-dimensional or three-dimensional free Bose gas, a Lorentzian 
momentum distribution being not normalizable in this case. 

To characterize the coherence properties of the gas we use correlation functions for the atomic field 
operator ^(x, y, z) in direct analogy with the correlation functions considered in optics for the photonic field 
& 

We define the first order correlation function as 

0i(*) = ($t(O,O,2)* (0,0,0)) (23) 

where the expectation value is taken in thermal equilibrium; this function is sensitive to the coherence of 
the atomic field between two points on the axis of the guide separated by a distance \z\. It can be written 
as the sum of the contributions of the transversally excited states and of the states in the transverse ground 
state. The two contributions behave in a very different way in the degenerate limit, that is in the limit 
H — ► 0. This is illustrated in figure § for a moderately degenerate regime. 

By calculations similar to the ones leading to Eqs.( |l4| , |l5| ) we find for the contribution of the transverse 
excited states to g\ in the thermodynamical limit: 



91 (*) = ^ \ 7/2 [ 1 - exp (- 2s ^ ±) " l \ 6X P ) • ( 24 ) 
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Figure 2: Correlation function gi(z) for the ideal Bose gas in the thermodynamic limit. The temperature is ksT = 

(c) (c) 

20hu>± and the linear density is py ln = l-lp lin where the critical density /9 lin is calculated with the approximate formula 
Eq.([l^). One finds numerically f3fi ~ —6.4 x 10~ 3 . Solid line: full correlation function, calculated numerically. Dashed 
line: contribution of the atoms in the transverse ground state of the trap, from the analytical formulas Eq.(|2l]) and 
Eqs.(|27|j2q). The coherence length is found to be l c ~ 3.5A. 



Using furthermore the fact that < Hujj_ <C ksT we set /x = and we expand the exponential function 
between square brackets to first order in P?ilu±: 



S>1 v 7 



(25) 



The maximal value of g^ is obtained in z = and is equal to p^ xis - The half-width of is on the order of 
0.75A. 

The contribution to g± of the states in the transverse ground state is given by 



(26) 



In the thermodynamical limit, we replace the sum by an integral. We use the Lorentzian approximation for 
the occupation numbers Eq.(]20|) and calculate its Fourier transform. This leads to a correlation function 
being an exponential function of \z\: 



(0) f . ™ji (o) z /i c 
91 [Z) - irh P ^ e 



with a coherence length 



lr 



* ,(0) 



(2m|/i|) 1 /2 2vr 



Pit 



(27) 



(28) 



The exponential decay of g\ found here is to be contrasted with g\ going to a finite value in the case of a 
three dimensional homogeneous Bose-Einstein condensate [[l?]]. The coherence length l c can be larger than 
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the thermal de Broglie wavelength A. It can be expressed in terms of the on- axis density using the relation 

Q (0) f Q) . (c) (0) . 

y 1 yz — U) — paxig — taxis' 

I ~ -^o (0) (29) 

A related issue is the expression connecting the on-axis and the linear density of the atoms in the transverse 
Bose-Einstein condensate: 

p£L „ 2C(5/2) 

>). - C(3/2) fiwjL >) ' W 

raxis I Im 

(c) 

This expression holds in the regime ksT 3> huj± so that a modest value of pn n above the critical value 
may actually correspond to a strongly degenerate regime p ax is^ 3 ^ 2.612. 
We define the second order correlation function of the atomic field as 



92{Z) 



{&(0, 0, z)*t (o, 0, 0)^(0, 0, 0)tf (0, 0, z)). (31) 



From a field point of view g 2 is the correlation function of the intensity of the field; it is a measure of 
the intensity fluctuations of the field. From a corpuscular point of view g 2 {z) is proportional to the pair 
distribution function of the atoms in the gas, that is the probability density to find a pair of atoms separated 
by a distance \z\ in the gas. 

For the ideal Bose gas in the grand canonical ensemble the use of Wick's theorem readily allows one to 
express g 2 in terms of g\ by performing all the possible binary contractions of the field operators: 

g 2 (z)=g 2 1 (0)+g 2 1 (z). (32) 

We find the unpleasant result that the gas is subject to large intensity fluctuations over a length scale on 
the order of the coherence length l c . In particular 

^ - 2 (33) 
alio)' 2 - (33) 

This value of two is typical of a bosonic bunching effect that manifests itself in an Hanbury-Brown and 
Twiss experiment in optics with thermal sources. It has to be contrasted with g 2 — g\ obtained with laser 
light for photons or with almost pure Bose-Einstein atomic condensates M. The output of our magnetic 
guide cannot be termed an 'atom-laser' if g 2 (z) significantly differs from g\ (0) over the coherence length of 
the field. Fortunately we shall see in the next sections that atomic interactions can improve the situation. 
Remark: 

A careful reader may argue that it is dangerous to use the grand canonical ensemble for the ideal Bose 
gas to calculate fluctuations of the number of particles, and therefore of the field intensity. This fear is 
justified when a condensate is formed, a well known problem in three-dimensions: the number of particles 
in the condensate has then unphysically large fluctuations. This problem does not take place here in the 
thermodynamical limit where no Bose-Einstein condensate is formed. More precisely one can deduce from 
g 2 the standard deviation of the total number of particles: 



A:V " V '" 2 (34) 



N \LJ 

The relative fluctuations in the number of particles become small as soon as the length L of the gas becomes 
much larger than the coherence length l c . 

To illustrate the strong intensity fluctuations of the ideal Bose gas in a dramatic way we introduce the 
Sudarshan-Glauber P representation of the many-body density operator [16|: 

a = Jv^\coh : ^)(coh : tf| P({*},{^*}). (35) 
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This expression is a functional integral J T>^> over the real part and the imaginary part of the c-number field 
^f(f ). It involves the coherent or Glauber state of the atomic field associated to 



|coh : = exp 



1 



d 3 r \^(f)\ 2 



exp 



d 3 f \I>(r )&( 



vacuum). 



Here the many-body density operator a is the grand canonical thermal density operator 

a oc exp[-(3(H - fiN)} 



(36) 



(37) 



where H is the Hamiltonian of the gas containing the kinetic energy and the trapping potential energy, and 
N is the operator total number of particles. The Glauber distribution function P can then be calculated 



exactly [16]. One expands the field ^ on the eigenmodes of the trap: 



LV2 



(38) 



where cj> n , n = 0, 1, . . . are the normalized eigenfunctions of the ID harmonic oscillator of frequency uj± and 
where the plane waves along z have a wavevector given by Eq.(||). Then the Glauber distribution function 
is simply a product over all modes of Gaussian distributions with squared widths given by the occupation 
number of the modes: 



P({tf },{**}) oc exp 



0:r 



(39) 



As this distribution P is positive the thermal equilibrium a can be viewed exactly as a statistical mixture 
of coherent states. One can then imagine that a given experimental realization of the Bose gas is a coherent 
state characterized by a field This field Vl/ is stochastic as it varies in an unpredictable way from one 
experimental realization to the other. As the coherent state |coh : v^) is an eigenstate of vE^r ) with the 
eigenvalue Vl/(r ) one can check that the correlation functions g\ and gi are equal to the following averages 
over all possible realizations of the c-number field: 



9l (z) = <**(*)*(0)> stoch 
g 2 (z) = <|*(z)| 2 |*(0)| 2 ) stoch . 



(40) 
(41) 



We have plotted in figure |3] the intensity of the field as function of position for two numerically generated 
realizations of Figure ||a corresponds to a non-degenerate situation; the only spatial scale for the intensity 
fluctuations is the thermal de Broglie wavelength A. Figure |3|b corresponds to a strongly degenerate regime; 
there are clearly two spatial scales for the intensity fluctuations, one on the order of A coming from the 
transversally non-condensed fraction and the other one on the order of l c due to the gas in the transverse 
ground state of the trap. The large intensity fluctuations at the scale of l c manifest themselves as "droplets" 
in the atomic density. 



3 The interacting case: a one-dimensional classical field model 
3.1 Model Hamiltonian for the transversally Bose condensed gas 

We have seen on the ideal Bose gas model that transverse Bose-Einstein condensation takes place in the 
magnetic guide at sufficient high density. Although we have not performed any detailed analysis we expect 
the same phenomenon to occur for the interacting Bose gas in the weakly interacting regime. 

At the presently considered low temperatures, interaction between the atoms takes place in s-wave 
mainly, and the relative wavevector of two colliding atoms is much smaller than the inverse of the scattering 
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Figure 3: Intensity of the atomic field as function of position along z axis for a given Monte Carlo sampling of the 
Glauber P distribution at thermal equilibrium. The temperature is kgT — 2QTilo^ and the size of the box is L = 400A. 
(a) Non transversally condensed regime /3fi = —0.1 corresponding to p\ in = 0.92p^. (b) Strongly degenerate regime 
fin = — 10~ 4 corresponding to p\ ln — 1.4pfc|. For (b) the coherence length is equal to l c /\ ~ 28. The critical density 
Piin * s calculated with the approximate formula Eq. ( |l3| ) . 



length a^d of the interaction potential. The interaction potential is then commonly replaced by an effective 
low energy interaction potential gsd&ifi 



r 2) 



19] with a coupling constant 



93d 



m 



-a?,d- 



The resulting Hamiltonian is expressed in terms of the field operator ^(x, y, z) as 



H 



3d 



dx / dy 



dz 



2m 



|grad + -g M ¥¥^ + (U (x, y) 



(42) 



(43) 



with the trapping potential given in Eq.(|l|) and with periodic boundary conditions along z. We have included 
for convenience the term —fiN where [i is the chemical potential so that the thermal equilibrium density 
operator is simply oc exp[— f3H 3( j\ m the grand canonical ensemble. 

In our magnetic guide geometry the situation is particularly simple if the typical interaction energy 
per particle p a xis<?3d is smaller than the quantum of transverse oscillation Huj±. Transverse Bose-Einstein 
condensation will then take place in a transverse wavefunction <p±(x,y) close to the ground state of the 
harmonic oscillator. If we wish to describe only the one-dimensional Bose gas of atoms in this transverse 
condensate, assuming that the remaining atoms have a much smaller density, we can neglect the contribution 
of the transverse modes to the field operator by setting 



V(x,y,z) ~ <j) ± (x,y)i/}{z). 



(44) 



This will eliminate the contribution of the atoms in the transverse excited states of the trap; we keep in 
mind however that these atoms are essential to ensure thermalization, with a temperature ksT > ?luj±. The 
reduced field operator if) obeys the usual bosonic commutation relations of a one-dimensional bosonic field. 
Inserting the approximate expression for ^ in the Hamiltonian results in the following model Hamiltonian 
for the transversally condensed Bose gas: 



H 



dz 



2m 



\dj\'- 



It corresponds to a one-dimensional Bose gas with a contact interaction potential g5{z\ 
particles with an effective coupling constant [ 20 ] 



9 = 93d 



dxdy\4>±{x,y)\ l 



(45) 

Z2) between 
(46) 
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3.2 Classical field approximation 

The model Hamiltonian Eq,( |45|) leads to an exactly solvable iV-body problem: exact eigenenergies and 
eigenfunctions are known [pl|, |22|. It is however not easy to extract information from the exact solution, 
even at zero temperature. 

We use here a simpler approach, valid in a sufficiently high temperature regime. The idea is to write 
the thermal equilibrium density operator (up to a normalization factor) as the result of a fictitious time 
evolution: 

!»— + (47) 

with the initial condition a(r = 0) equal to the identity operator. This evolution corresponds to the so-called 
imaginary time evolution. It leads to a density operator at 'time' r given by: 

a = e' T& (48) 

which is (up to a normalization factor) the thermal equilibrium density operator at temperature /cgT = 1/r. 
High temperatures correspond to low values of r that is to short 'time' evolution. 

To take advantage of the fact that the 'time' evolution is short it is better to rewrite Eq.(|47|) using some 
of the representations of the density operator introduced in quantum optics. We use here the Glauber P 
distribution already introduced in section |2|. It is defined as in Eq,([35|) with the difference that ip is now 
a function of the coordinate z only. So we apply the Glauber transform to both sides of Eq.(^Tj). The 
transform of the right-hand is performed with the following rules: 

Glaub [i>(z)v] = V(^)Glaub [i>] (49) 
Glaub[^ f (z)z)] = (ip*(z) Glaub [u] (50) 

where v is any operator and Glaub [v] stands for the Glauber P distribution of v. The first rule Eq,(f4"g|) 
comes from the fact that the coherent state |coh : ip) is an eigenstate of ip(z) with the eigenvalue ip(z). The 
second rule Eq.(|50|) involves a functional derivative with respect to the field ip, the fields tp and ip* being 
formally considered as independent variables. Its derivation closely follows the one for single mode fields in 
One first uses the following identity: 

d^ z) [|coh : ^)(coh : ip\] = (ft(z) - ip*(z))\coh : ^)(coh : ip\ (51) 



as can be checked from the definition Eq.(|3q) transposed to the one-dimensional case. Then one integrates 
by part in the functional integral over the field to convert the derivative over the dyadic |coh : ip){coh : ip\ 
into a derivative of the Glauber P distribution. We finally obtain the Fokker-Planck type equation for the 
fictitious time evolution of the Glauber distribution: 

d T P = -EP- J dz [^(F(z)P) + |^(^ 2 P)+c.c] . (52) 



The first term in Eq.(|52|) involves a multiplication of P by a functional E which is simply the Gross- 
Pitaevskii energy functional |17| , obtained by replacing in the Hamiltonian the field operator by the c-number 
field tp: 



E[{iP},{iP*}]= [ dz 
Jo 



(53) 



|-I^| 2 + ^| 4 -^P 
2m 2 



If this term was alone in the evolution equation for P we could readily integrate it to obtain: 

Pci ass [m,{r}]=e- TE ^ r ^. (54) 

We have termed this solution -P c i ass as it corresponds to the thermal Boltzmann distribution for a classical 
field ip with an energy given by ^[{'0}, {V'*}]' This is nineteenth century equilibrium physics for fields. 
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The next term in Eq.(p2[) can be termed a force term by analogy with the Fokker-Planck equation, as it 
involves a first order derivative in ip. The 'force' functional is given by: 

-,2 

*{z). (55) 



2 f (*) LLij' l ~ jj 2 



2m 



At sufficiently short 'time' r, that is at high temperature, the effect of the force term during the evolution 
'time' r is to shift the field ip(z) from its initial (random) value by the amount F(z)t. One may hope that 
this shift is negligible for high enough temperature. 

Let us calculate this shift for the ideal Bose gas. The field ip is then conveniently expanded on plane 
waves with momenta given by Eq.(|^): 

^( z ) = z^2 ak *JT/2- ( 56 ) 

k z 

Taking now k z and a kz as coordinate and field variables (rather than z and ip(z)) we can write: 

E = J2 £ k^* kz a kz and F(k z ) = --d a *^E = --e kz a kz (57) 

k z 

where we have introduced the mode eigenenergy e kz = h 2 k 2 J{2m) - /x. The shift F(k z )r has a negligible 
effect on field mode k z if it is small as compared to a kz ; this leads to the condition 

k B T^e kz . (58) 

The classical field approximation -P c i ass is therefore an acceptable approximation for the modes with an 
energy much smaller than /cgT. The occupation number of such modes can then be obtained from the 
classical energy equipartition formula: a mode of a complex field corresponds to a one-dimensional harmonic 
oscillator so it has a mean energy equal to fc^T, and 

(\a k f) = ^. (59) 

e k z 



This formula coincides indeed with the quantum Bose formula Eq.(y) in the limit Eq. (|5q) , that is in the 
limit of a large occupation number of the mode. 

So when can we use the classical field approximation Eq.(|54|)? The answer depends on the observable 
quantity we wish to calculate. 

For the calculation of the mean energy the classical field approximation is never acceptable: in the 
absence of energy cut it predicts an infinite mean energy, the well known blackbody radiation catastrophe. 
To save the situation one has to introduce an energy cut e cut on the order of kgT to reproduce 'by hand' 
the fact that the Bose formula gives an exponentially small occupation number to modes with eigenenergy 
much larger than k^T . The mean energy is then finite but depends on the precise value of e cu t- 

The conclusion is different for the calculation of the correlation functions g\ and g2- We have actually 
already used the classical field approximation in the derivation of gi, see the approximation Eq.(pO|)! This 
did not lead to any divergence, a fortunate feature peculiar to the free one-dimensional Bose gas. The 
classical field predictions for g\ and 52 therefore do not depend on the energy cut e cut provided that the 
energy cut is large enough. For the ideal Bose gas the condition is e cu t 3> as the kinetic energy width 



of the Lorentzian Eq.(20) is as e cut ~ fceT the classical field calculation of g\^ requires ksT S> we 
recover Eq. fllif ). 



The validity conditions of the classical field approximation for the interacting case are more subtle to 



derive and will be discussed in §3.5. 



Finally, the last term in Eq.(52) can be termed a diffusion term by analogy with the Fokker-Planck 
equation, as it involves second order derivatives in ip. The corresponding 'diffusion' matrix in point z, given 
by: 

^-?(.l (60) 



4 V V 
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is however a non positive matrix: one can check that the field quadrature along ip is squeezed by the 'diffusion' 
while the field quadrature orthogonal to ip gets anti-squeezed. This non-positivity of the 'diffusion' matrix 
makes it impossible to perform a stochastic, Brownian type simulation of Eq.(|52]), which would have provided 
an exact numerical solution to the problem. 

3.3 How to calculate the correlation functions in the classical field approximation 

We take here as Glauber P distribution the classical approximation Eq. (|54|) , without introducing any energy 
cut, and we wish to calculate the stochastic averages Eqs.( pC|j4l] ). This amounts to calculating ratios of 
functional integrals over paths parametrized by z £ [0, L].. For example the first order correlation function 
of the field ^f(z) in the approximation Eq.(E3) is given by 



9l (-)HMo.o)|* Jd ^^ viie . mithti! , n ■ 

In the above formula the functional integrals are performed over all possible closed paths, as the gas is 
subject to spatially periodic boundary conditions; we have split the functional integral as a regular integral 
over the value of the path tp(0) in z = and a functional integral over the set Tmq) of all paths starting 
with the value tp(0) in z = and ending with the same value in z = L. We explain here how to calculate 



these functional integrals. The reader not interested in technicalities may jump directly to §3.4. 

To calculate functional integrals like Eq.(|6l"|) it is of course possible to use a Monte Carlo method to 
sample the distribution of ip, e.g. by representing on a finite spatial grid with step dz and by evolving ip 
one time step dt after the other according to the stochastic evolution 



dip(z) = —dt 



2m 



Hz) + (^y\(z) (62) 



where the d£(z)'s are statistically independent complex noises of variance 2dt. 

There exists however a more elegant and much faster solution. One can use the link between path 
integrals and quantum mechanics propagator put forward by Feynman. The functional integral over the 
classical complex field then corresponds to a propagator in imaginary time of the quantum mechanical 
problem of a particle in two dimensions. The point by point analogy between the two problems is specified 
in the translation table: 

classical field problem quantum mechanical analogy 

path integral *-* quantum propagator 

abscissa z <-> time t 

Re (i/j(z)) <-> position x(t) 

Im(^>(z)) *-* position y(t) 

V* e -/"WM*-}] _> ( Xf , y f \e-*W*\xi t Vi ) . 

We have to identify the Hamiltonian H of the equivalent quantum mechanics problem. We postulate 
the following form: 

K = P -^ + V(x,y). (63) 



The imaginary time propagator is then expressed in terms of the path integral [23]: 



(x f ,y f \e- tn / h \x l ,y t )= / Vx(r) Vy{r) e - s K*M«/>]/ fi (64) 

Jx{o)=Xi Jy(0)=yi 
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where the action S is a functional of the path x(r),y(r) 



S[{x},{y}}= fdr 
Jo 



(65) 



One identifies this action with %(5E and one uses the translation table to obtain the values of the parameters 
of the equivalent quantum mechanics problem: a mass 

u - 1st (66> 

and a potential 

V(x, y) = h(3 [|(x 2 + y 2 ) 2 - ^{x 2 + y 2 )] . (67) 

This potential is rotationally invariant, as a consequence of the U(l) symmetry of the Hamiltonian Eq.(|45|). 
We can therefore classify the eigenstates of V with two quantum numbers, an angular momentum m and a 
radial quantum number n = 0,1, . . . . We call 4>™ the corresponding normalized eigenvector with eigenvalue 
e™. As usual the absolute ground state of H is of angular momentum m = and radial quantum number 
n = 0. 

We translate the functional integrals of Eq.(^) into quantum propagators. In particular we note that the 
integral over ip(0) in Eq.(^TJ) corresponds to an integral over all possible initial coordinates of the particle, 
that is to a trace over all possible initial quantum states of the particle. We finally obtain: 



9i(z) = |0± (0,0)| 



Tr [e-W/h] 

and a similar expression for g 2 . 

Physically, as there is no Bose-Einstein condensation along z, one expects that the length of magnetic 
guide L in the experiment is much larger than any correlation length of the gas. One can then take 
the thermodynamical limit along z, putting L to infinity while keeping a constant chemical potential fi 
(this ensures that the mean linear density N/L tends to constant) |24|]. In this case exp[— LTi/fi] becomes 
proportional to the projector on the absolute ground state of 7i: 

e -LH/H „ e _ Le ^o /?w=0) w=0| _ (g9) 

The thermodynamical limit approximation greatly simplifies the expressions for the correlation functions, 
as the trace Tr can be restricted to the ground state of H: 

9i{z) = |^(0,0)| 2 (C =0 |(^-%)e- 2(W - e ^° )/?i (x + %)|C =0 ) (70) 

g 2 (z) = |^(o,o)| 4 (C =0 l(^ + y 2 )e" 2(H ~ £ ^° )/?i (^ + y 2 )IC =0 )- (71) 

The operator x + iy maps the absolute ground state to a state with angular momentum equal to unity. If 
we restrict for simplicity to the large z limit, the operator e~ z ^ in Eq,(f70|) becomes proportional to the 
projector on the ground state <^(j 1=1 of TL with angular momentum m = 1 so that 

g x {z) ~ |^(0,0)| 2 a 1 exp(-K 1 z) (72) 

with 

K 1 = ^r 1 -£^ =0 )/n and a 1 = \{^= 1 \x + iy\ ( j^= )\ 2 . (73) 

A similar analysis can be applied to the correlation function g2, with the difference that the operator x 2 + y 2 
in Eq.([n|) maps the absolute ground state to a rotationally invariant state. In the large z limit we keep the 
contributions of the first two eigenstates with m = to obtain 

g 2 {z) ~ 5l 2 (O) + |0 ± (O,O)| 4 a 2 exp(-/ t2 z) (74) 

with 

K2 = (er°-eF*)/h and a 2 = \{^=°\x 2 + y 2 \^=®)\ 2 . (75) 
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3.4 Results of the classical field approximation 



We wish to calculate the correlation functions g\ and 52 using the formalism of § |3.3| . One has then to solve 
the quantum mechanics equivalent problem of a particle in two-dimension with the Hamiltonian Eq.(|63|). 
This Hamiltonian can be diagonalized numerically. We wish to express the results in terms of the linear 
density of the transversally condensed Bose gas, rather than in terms of the chemical potential [i. For a 
given linear density we therefore have to adjust p in order to satisfy 

pfl = <^(0)&0)>. (76) 

The problem can be simplified by an efficient parameterization. We express the coordinates x and y in the 
quantum mechanics analogy (which correspond to the real and imaginary part of ip) in units of the square 
root of the linear density. We express the physical length z in units of p$\ 2 /(2tt) where A is the thermal 



de Broglie wavelength Eq.(lO). We then find that once p has been eliminated there is a single parameter 



left in the classical field theory: 

x _&if^l (77) 

X k B T 2vr ' { 1 

We plot in figure || as function of \ the coefficients «i 2 giving the long distance behavior of the correlation 
functions as defined in Eqs.(73,75|). In the limit of a vanishing x we recover the results of the ideal Bose 



gas, Eq. (P8[) , with k\ = l/l c , and Eq. (|32|) leading to K2 = 2k\. For an increasing interaction strength, x 
increases: the coherence length 1 / k\ has a modest increase by up to a factor two; the correlation length of 
the intensity fluctuations l//«2 is dramatically reduced by the atomic interactions and becomes much smaller 
than the coherence length, a positive point already ! 

What happens to the amplitude of the intensity fluctuations ? In an ideal 'atom-laser' there is no 
fluctuation of the field intensity; the deviation from this ideal situation can be measured by the ratio of the 
maximal to the minimal value of g2, that is 

92(0) = (4 m=0) |(^ 2 + ^ 2 ) 2 l4 m=0) ) m 

A® (<4 m =V + ^ m = 0) >) 2 ' 

We have plotted this quantity in figure |j| It is equal to two for the ideal Bose gas, as predicted by Eq.(^). 
It is sharply reduced by the atomic interactions for low values of x then it slowly goes to unity for large 
values of x- 

Can we understand the origin of the reduction of intensity fluctuations using the quantum mechanics 
analogy ? According to Eq.([Tq) this amounts to understanding the fluctuations of the distance of the 
quantum mechanical particle from the origin in the ground state of 7i. By inspection of Eq.(|67|) giving 
the trapping potential seen by the quantum mechanical particle we realize that there are two situations 
depending on the sign of the chemical potential. For a negative chemical potential, as in the case of the 
ideal Bose gas, the potential V has an absolute minimum in x = y = so that |<^™ =0 | 2 is localized around 
the origin (see figure |6|a): large fluctuations of the intensity of the field are expected. For a positive chemical 
potential, which is the case for strong enough repulsive interactions, the potential V(x, y) has a Mexican 
hat shape: it has a local maximum at the origin and a global minimum on a finite circle (see figure ||b). 
The ground state wavefunction then tends to be localized around the circle. The critical regime for the 
apparition of the Mexican hat potential is such that p = 0; we find numerically that this corresponds to 

Xc * 0.28. (79) 

This low value of x explains why a sharp variation is obtained at the scale of figure ||. 

In the large x limit the quantum mechanical particle will get more deeply bound at the bottom of 
the Mexican hat potential. Writing Schrodinger's equation for <ft™~ 0,1 in polar coordinates and treating 
perturbatively the deviation of the Mexican hat plus centrifugal potential from a harmonic approximation 
we obtain after some algebra the large x's expansions: 
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We note that the correlation length of the intensity fluctuations l/«2 becomes proportional in the large 
X limit to the so-called healing length £ of the gas, a crucial parameter in the theory of Bose-Einstein 



condensates [17 



1/2 



(0) 



y/2' 



(81) 



3.5 Validity conditions of the classical field approximation 

The validity condition of the classical field approximation in general depends on the observable to be 
calculated. Here the relevant quantities are the first and second order correlation functions gxp- 

Let us recall briefly what happens in the ideal Bose gas case. The correlations functions are then 
characterized by a single length, the coherence length l c = 1/kx = 2/«2. The plane wave modes of the 
field contributing to gx,2 have therefore wave vectors on the order of Kx or less; for these modes the validity 
condition Eq.fl58j) reads &bT 3> e Kl . The mode eigenenergy e Kl is equal to Ti 2 K\/(2m) — \i\ using the value 
Eq.(|28|) of the coherence length l c = 1/kx we arrive at 

fc B T»|/z|. (82) 



This result was already obtained in Eq.(lC). It can be rewritten with Eq.(21) as a hierarchy among the 
three relevant lengths of the problem: 



fC) <A<l c = «i 1 . (83) 

What happens in the interacting regime x > 1 ? A difficulty is that the Hamiltonian for the classical 
field ifi in Eq.(|53|) is no longer quadratic in ip so that it is not straightforward to calculate energy eigenmodes. 
Can we find some good quadratic approximation to it? 



One could think to use the Bogolubov approach [17]. In this approach one identifies the field ipo min- 
imizing the energy Eq.(|53"D; in the present homogeneous case ipo is z-independent, ipo = (fi/g) 1 / 2 exp(i9o) 
where 9$ is an arbitrary constant phase. Then one splits the field as vp(z) = tpQ + 6ip(z). Under the assump- 
tion of \6ip\ <C |^o I one neglects in Eq.(53) the terms cubic and quartic in Sip which leads to a quadratic 



Hamiltonian that can be diagonalized. This approach is not well suited to the present situation in the large 
L limit; it predicts 

W\ 2 e" lL - (84) 

Although the emergence of a coefficient proportional to n\ is promising, the Bogolubov quadratization 
procedure is not justified when L exceeds the coherence length of the gas. This is physically not surprising: 
the finite range of first order field coherence in our one-dimensional geometry is precisely due to large 
fluctuations of ip(z) away from i/jq. 

Fortunately it is possible to adapt Bogolubov's idea taking advantage of the weak intensity fluctuations 
of the field in the large x regime. We split the field in a modulus and a phase factor 

xj}{z) = p^^zy 6 ^. (85) 
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Figure 4: In the one-dimensional classical field model, the inverse field coherence length k% and the inverse intensity 
correlation length K2 as function of the universal parameter \ defined in Eq.(|7?|). 
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Figure 5: In the one-dimensional classical field model, indicator of the intensity fluctuations of the field as function 
of the universal parameter \ defined in Eq. ( |77| ) . 
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Figure 6: In the quantum mechanics equivalent to the classical field problem, potential V(x, y) seen by the quantum 
mechanical particle, (a) For a low value of \ = 0.1 and (b) for a large value of \ = 20. The units on the axes are such 
that V depends only on the parameter \- 



We recall that the intensity p(z) and the phase 9(z) of a field are Hamiltonian conjugate variables in the 
same way as ip(z) and ip*{z) are. The intensity p{z) has only small fluctuations away from the most probable 
value po = \ip \ 2 = p/g: 



p(z) = po + 5p(z) with \Sp(z)\ <C po- 



(86) 



The mean value of p(z) is simply the linear density of atoms in the transverse ground state p{^; to lowest 
order, po can be identified with this mean density. The phase 9 has on the contrary large fluctuations 
away from any fixed constant phase 9q over distances larger than a few coherence lengths. So contrarily to 
Bogolubov's method we do not assume 9{z) — 9q to be small. 

The potential energy density of the field g\tp\ 4: /2 — p\ip\ 2 leads to a constant term and a term quadratic 
in Sp. In the kinetic energy density of the field we replace p(z) by its lowest order approximation giving 
a non-zero contribution. Up to a constant, this leads to a quadratic approximation for the Hamiltonian 
Eq.( 
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(87) 



One has finally to diagonalize this quadratic Hamiltonian using canonically conjugate variables. An easy 
way is to expand the field variables on the eigenmodes of the linear time evolution equations: 

-■2 



fid t 5p{z) 
-hd t 9(z) 



quad 



m 



p d 2 z 9(z) 



ds P (z)Equa.d 



Ti 2 



4m po 



d 2 z 5p(z) +gdp(z). 



(89) 



These well known equations have the form of linearized hydrodynamics equations for a superfluid 
Eigenmodes are plane waves of wavevector k z obeying Eq.(||) and with eigenenergy 



£ k z 



h 2 k 2 (h 2 k 2 



2m V 2m 



+ 2po5 



1/2 



(90) 



This is the famous Bogolubov spectrum; the field energy will then formally appear as a sum of decoupled 



harmonic oscillators, corresponding to an ideal Bose gas of quasi-particles 17]. 

We can finally reproduce the reasoning performed above in the ideal Bose gas case to identify the 
validity condition of the classical field approximation. The temperature must be much larger than the mode 
eigenenergies at wavevectors k z = k%,k,2 relevant for the correlation functions g\,g2- 



k B T^> 



(91) 
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As K2 3> K\ in the large x regime, the condition involving K2 is the most stringent one. Prom Eq.(fl]) we 
find formally the same condition as the ideal Bose gas case 



k B T > fi (92) 
with now a different expression for the chemical potential, [i ~ P\^9- We note that in the large x regime 



and the high temperature regime Eq.(92) we have the following hierarchy among the various relevant scales 
of the problem: 

(pnl) 1 < A < «2 1 - £ < lc = < « (93) 



where the 'one-dimensional scattering length' a is defined as 

g = —. (94) 



ma 



The property that the microscopic scale a is larger than the other physical lengths of the problem (with the 
exception of L of course!) ensures that the one-dimensional Bose gas is a weakly interacting Bose gas. The 
fact that this is a necessary condition to a classical field approximation (like the Gross-Pitaevskii equation) 
is known for three-dimensional Bose gases, with the difference that the three-dimensional scattering length 
a^d has then to be smaller than the macroscopic scales of the gas p7| . In particular the small gaseous 
parameter is a^/^, for three-dimensional Bose gases whereas it is £/a for one-dimensional Bose gases J21[ | . 
This is not surprising if one realizes that the weak interaction limit g —* leads to a — > 00 in one-dimension. 



3.6 Comparison with results in the literature 

The model of the one-dimensional interacting Bose gas has already been studied by several authors without 
performing the classical field approximation. 

A first line of thought deals with the exact solvable model of a contact interaction potential [^l], p2| . 
To our knowledge a full calculation of gi(z) and #2 ( z ) a t finite temperature has not been performed. The 



value of (72(0) can be found in [21| at zero temperature. Exact results have been obtained in the strongly 
interacting regime g — > +00: the correlation function gi(z) has been calculated in at zero temperature; 
the correlation functions gi(z) and <?2 (z) at finite temperature have been obtained in |27[| . These regimes 
are different from the finite temperature, weakly interacting Bose gas considered in the present paper. The 
extension of the calculation of 52(0) to finite temperature seams feasible but we do not know any reference. 

A second line of thought is to consider the regime of low intensity fluctuations of the field (our large 
X regime): one may take advantage of the weakness of the intensity fluctuations by a linearization of the 



equations of motion in the hydrodynamic point of view |2q, 29, 30], a Bogolubov type approach [31], a 
quadratization of action in a path integral formulation [p^, |32[] . 

We have checked that the classical field predictions in the asymptotic limit x>l reproduce the results 
of e.g. [31] when (/jbT) 1 / 2 3> {p^g) 1 ^ 2 ■ The advantage of the classical field approximation is that it is not 



restricted to the large x regime so that the transition from the ideal Bose gas to the strongly interacting 
case can be studied. 

In preparing this article we have discovered that a one-dimensional classical field model very similar to 
our model has been studied in with a different physical motivation. 



4 Conclusion and perspectives 

In this paper we have discussed several aspects of our proposal for the production of a continuous 'atom-laser' 
source, consisting in evaporatively cooling an atomic beam in a long magnetic guide. From the classical 
Boltzmann equation and for expected typical parameters we have estimated the length required to reach 
the quantum degenerate regime, a few meters for two-dimensional evaporation; the corresponding loss on 
the atomic flux is only two orders of magnitude. 
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We have also characterized the coherence properties of the output beam once quantum degeneracy is 
reached. The gas is expected to experience transverse Bose-Einstein condensation in the guide, leading to a 
transversally monomode output beam. We have therefore introduced a one-dimensional model for the gas, 
that we have solved in a classical field approximation. The coherence length of the field along the axis of 
the magnetic guide can be much larger than the thermal de Broglie wavelength, typically by one order of 
magnitude. The intensity fluctuations of the beam are very large for an ideal Bose gas but are strongly 
reduced by repulsive atomic interactions when the healing length £ of the gas becomes smaller than the 
coherence length of the gas. 

Possible extensions of this work are the discussion of superfluidity properties of the 'atom-laser' and a 
complete three-dimensional modeling of the interacting gas in the magnetic guide. 

We acknowledge useful discussions with Gora Shlyapnikov, Philippe Grangier, Gordon Baym and Tony 
Leggett. We thank Alice Sinatra for helpful comments on the manuscript. Laboratoire Kastler Brossel is a 
unite de recherche de l'Ecole normale superieure et de l'Universite Pierre et Marie Curie, associee au CNRS. 
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